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ABSTRACT 

Time evolution of tunneling phenomena in medium is studied using a standard 
model of environment interaction. A semiclassical formula valid at low, but finite 
temperatures is derived in the form of integral transform for the reduced Wigner 
function, and the tunneling probability in thermal medium is calculated for a gen- 
eral tunneling potential of one dimensional system. Effect of dissipation, its time 
evolution in particular, depends on the behavior of the potential far beyond the 
barrier. 
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Tunneling phenomena are genuinely a quantum effect. When they occur in some 
surrounding medium, the important question arises as to whether the tunneling rate 
is enhanced or suppressed by the environment effect. There are already many works 
on this subject and most past works deal with a system in equilibrium as a 
whole. The Euclidean technique such as the bounce solution is often used in this 
context [0,0]. Our approach here is different, and we attempt to clarify dynamics 
of time evolution starting from an arbitrary initial state of the tunneling system, 
which can be either a pure or a mixed state. Moreover we find it more illuminating 
to use a real-time formalism instead of the Euclidean method much employed in the 
literature p. 

Although one can think of many applications of our general framework here, ap- 
plication to cosmology is our prime target 0. The cosmological environment is in 
many aspects the simplest since the environment is almost spatially homogeneous 
and for many purposes its state can be characterized by a single temperature pa- 
rameter. Furthermore cosmic evolution gives rise to decrease of the temperature, 
dropping some degree of freedom out of equilibrium with the rest. This gives a 
natural setting for separation of a subsystem from the environment, which is at the 
heart of our approach. 

The standard model of environment and its interaction with a one di- 

mensional system which we describe by a potential V{q) introduces an infinitely 
many, continuous harmonic oscillators of environment and a bilinear interaction, 
q J duo c{uj)Q{uj) . Here Q{uj) is the environment oscillator coordinate of frequency uo 
and c{(jj) gives a coupling strength of the system-environment interaction. Although 
the form of interaction seems rather limited, it is generally believed that details of 
the environment interaction with the system should not be important and this form 
should suffice. With the total system specified, dynamics is given by the quantum 
equation of motion. 

Quantum Langevin equation is derived ^ by eliminating the environment vari- 
able Q{u! ,t); 

5 + ^ + 2 dsajit - s)qis) = Fg(t) , (2) 

where Fgit) has an explicit, linear dependence on the initial values, Qi^uj) , Pi{uj) = 
Qi{uj). In a thermal bath of environment temperature T = l/j3 the random force 
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from the environment is characterized by the correlation function, 

({^q(^) >^Q(-5)}+)cnv = / rfu;r(u;) cosu;(r - s) coth(— ) , (3) 

with r{uj) = {uj) / {2uj) . The kernel function «/ in eq.(0) is related to the real-time 
thermal Green's function and is given by a/(t) = — dujr{ijj) sin(co't) . 

An often used simplification is the local, Ohmic approximation taking r{ijj) = 
rjLj/n with Uc = 0, which amounts to a;/(r) = 6uj'^6{t) + ?7 5'(r) . This gives 6uj^q + 
T]q in the Langevin equation (0). The parameter 600"^ is interpreted as a potential 
renormalization effect or the mass renormalization in the field theory analogy, while 
rj is the Ohmic friction coefficient. This approximation breaks down at early times 
IQ, but it is useful in many other cases. 

Our discussion starts from the master equation for this system, which is written 
for convenience in terms of the Wigner function fwiq ,p ,Q{i^) , P{i^)), a Fourier 
transformation of the density matrix p{q , q' , Q{lj) , Q'{uj)) with respect to the relative 
coordinate, q — q' , Qiyj) — Q'{uj), 

dfw dfw f . ( j^, ^ dfw , ... d , ^/ x ^ w \ 



dt dq J \ dQ{uj) dP{uj) dp 

The crucial observation for subsequent development is how the Planck constant 
h enters in the master equation @). In one dimensional quantum mechanics the 
semiclassical approximation is excellent when the potential barrier is large, and we 
assume that this is also true in the presence of the system-environment interaction. 
In the semiclassical — limit we have 

myy^^d-pj-^y-^d-pji^'^^ii^^- 

The resulting equation, being identical to the classical Liouville equation, has an 
obvious solution; 

fw{q,P,Q ,P) = J dqidpi j dQidPi fl^ (g^ ,Pi,Qi, Pi) 

■S {q - gel) S {p - Pel) S {Q - Qci) 6 {P - P^) , (6) 

where qci{qi ,Pi , Qi , Pi;t) etc. are the solution of (HD, taken as the classical equation. 

We consider the circumstance under which the tunneling system is initially in a 
state uncorrelated to the rest of environment. Thus we take an un correlated initial 
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state of the form 

p(i) = ^ pW ^ 

to get the reduced Wigner function after the 

trivial Q{uj) ,P{uj) integration, 

K{q ,p,qi,pi;t) = J dQidPi fw,QiQi , Pi) Hi - Qd) Hp - Pa) ■ (7) 

The problem of great interest is how further one can simplify the kernel function K 
here. 

In many situations one is interested in the tunneling probability when the envi- 
ronment temperature is low enough. At low temperatures of T ^ a typical frequency 



or curvature scale Ug of the potential V, one has Ug y Qi{i^) , y P^{^) = C'[-\/T] <^ 
^yuJ^ . Expansion in terms of Qi{u!) , Pi^uj) is then justified. Thus, we use 



Hq-Qci) = f — ^ exp [ z Ag ( g - gel ) ] 
J Ztx 



dXa 

— - exp 



valid to the first order of Qi{uj) ,Pi{uj). A similar expansion for 6 {p — Pc\) using 
p^^\ p^^\uj) , p[i\uj) , also holds. 

Gaussian integral for the variables Qi{uj) , Pi{uj) ,Xq,Xp can be done explicitly 
with eq.(|^), when one takes the thermal density matrix for the initial environment, 
Pq' . The result of this Gaussian integral leads to an integral transform of the Wigner 
function, f^^ —>■ using the kernel function of 



K{q,p,qi,pi;t) = — exp 



1/ (0) (o)x ^ ( Q -Qd 

7;{<1 - Qci ^ P - Pel ) J (0) 

^ VP -Pel 



(9) 



where the matrix elements of JT" ^ = lij are given by 

1 ^ (3uj 1 ^,,2 



hi = - du coi\i^-\z{uo,t)\\ (10) 

1 Qlo 1 

/22 = - / c/u;coth^-|i(cu,t)p, (11) 

r doj cot\i^-^{z{to ,t)z*{uj ,t)) , (12) 
z Juic z to 

where z{uj ,t) = q''^\uj ,t) + ion q^i\uj ,t) , and z{uj ,t) = p^^\uj ,t) +iuj p''^p{uj ,t) . 

Quantities that appear here ought to be determined by solving differential equa- 
tions with a given initial condition; the homogeneous Langevin equation and its 
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linearized equation of the form 




+ 2 [ ds ai{t - s)q^\s) =0 



(13) 



dh{u,t) (d^V\ 



z{(jj,t) + 2 / ds aiit — s)z{(jj , s) 



c{uj)e' 



(14) 



dt'^ \ dx"^ j (0) 







The same hnear equation as for z{uj , t) holds for z{uj , t) . The initial condition is 



q(^>(t = 0) = Qi, p^^> (t = 0) = Pi , z{uj ,t = 0) = , z{u ,t = 0) = . Once 



the reduced Wigner function is known, one can compute any physical quantity by a 
phase space integration. 

The physical picture underlying the formula for integral transform should be ev- 
ident; the probabihty at a phase space point (g ,p) is dominated by the semiclassical 
trajectory (environment effect of dissipation being included for its determination) 
reaching {q,p) from an initial point {qi,Pi) whose contribution is weighed by the 
quantum mechanical probability f^) initially given. The broadning of the trajectory 
due to the environment interaction is given by the width factor J hj- 



To proceed further, we distinguish two cases of the potential type, depending on 
the value of V{oo) relative to the local minimum value Vm in the potential well. The 
kind of potential we have in mind is depicted in Fig.l. When V{oo) < Vm (Fig. la), 
the classical motion is usually monotonic ending at x = oo, while for V(oo) > Vm 
(Fig. lb) the motion is a damped oscillation towards Xq at the true minimum, unless 
the friction is large. If the friction is larger than a critical value of ~ 2u!q with 
ujq given by the curvature of the potential at the global minimum, there occurs 
the overdamping such that q'^^ —>■ Xq monotonically. A typical interesting case for 
V{oo) > Vm is the asymmetric double well as may occur in the first order electroweak 
phase transition [Q. 

In the rest of this paper, we consider a few problems to illustrate consequences of 
our general formula of the integral transform. One problem is the outgoing flux to 
the overbarrier region of the type of potential of Fig. la, assuming an initial energy 
eigenstate under the whole potential V. The other problem is the tunneling proba- 
bility for the type of potential of Fig. lb. We take here an initial thermal state of the 
same temperature as the environment, given by the density matrix oc X)n e~^'^" 1^) (^1 
with |?T,) the exact energy eigenstate under V. This choice is made in order to mimic 
the cosmological first order phase transition in which the local minimum in the left 
was the true minimum above a critical temperature. We insist on the energy eigen- 
state because it is the most stable state, being stationary when the environment 
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interaction is switched off. In this sense the choice of the energy eigenstate and its 
mixture is useful in many apphcations. 

We first discuss the tunneling probability for the case of V{oo) < Vm- Suppose 
that we are primarily interested in the tunneling probability from a state localized in 
the left well into the right region far beyond the barrier. For discussion from further 
on we assume uj^ ^ V/j, where cu* is the frequency defined at the potential bottom 
in the well and Vh the barrier height seen from the bottom. Since higher energy 
states contribute roughly with the weight factor, ^-'^^'yH-^'Pi)!'^* ^ where , "pi) is 
the hamiltonian of the hypothetical oscillator of frequency u;* without the barrier 
in the right, the high momentum component above the barrier is suppressed by 
g-2yft/w*_ assume that this factor is much smaller than a typical barrier pene- 
tration factor \T{E)\^ in the overbarrier region, usually given by the WKB formula, 
exp[— 2 / dx\p{x) \ ] . In this way we may restrict the phase space region in the inte- 
gral transform to the region D defined by |Vo| < -|- |a;o (g^ — xq)^ < -|- |Vo| , 
with Vo the minimum potential value at Xq. 

The flux at x is calculated from an integral of the form, ~ / dpp f^\x ,p;t) , 
giving 



f 1 

I{x,t) = dqi dpi {qi,Pi) ■^y^f = exp 



{x-qS^r 



21 
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Pr + g(--aj. (15) 

The initial Wigner function fl^{qi ,Pi) is a Fourier transform of ip*{qi — 2)'^{Qi + 1) 
with respect to ^, where ip is the initial wave function. The phase of the wave 
function plays an important role for the region, qi > where x^, is a turning point 
in the overbarrier region, and one may use the stationary phase approximation for 
the ^ and Pi integration. We find a stationary point at ,^ = , Pi = I{q-i) , = 
^ {'ip*{x)ip'{x) — ip' *{x)ip{x)) /\ip{x)\^ . The quantity I{x) is the usual quantum 
mechanical flux associated with the initial quantum state. 

In the region of x ^ there is always a classical trajectory that reaches the point 
X of the Gaussian peak from an initial qi in the range q^ > x^,, and the entire region 
within the Gaussian width y/Tn is fully covered by the integral. The factor outside 
the Gaussian exponent is well approximated by its peak value in the weak coupling 
case. Thus, one obtains I{x,t) ~ |'0(xc)|^ (~^c) , where one determines Xc{x ,t) 
from qS\xc , I(xc) ,0,0;t)^x. Using the WKB wave function for the energy eigen- 
state at qi > x„ i){qi) = T(£;) exp [i p(x)] /^/piqi), P(x) = \/'^(E-V(x)) , 
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one has a factorized form of the flux; 



/(x,t)^|T(E)|V(x,t;E), f = '^^[^^^^^ ■ (16) 

We note that the potential renormalization effect should be included for |r(£')p, as 
emphasized in 0. 

This formula for the flux reproduces our previous result for a specific potential 
of the inverted harmonic oscillator [|TU|, V{x) = —^uj^x^. In this case the Gaus- 



sian integral gives the exact result, relying neither on the low temperature nor on 
the stationary phase approximation. The present approach actually improves our 
previous result; 

^ ^ 9 + gljxc) _^ 9 +ujBg 
gl{xc) + guj^Xc ujB{g + ujBg) ' 

where tu^ ~ — ri/2 is a diagonalized frequency, and the initial flux 



I{xc) = \J2{E — V{xc)) in the WKB approximation. The function g{t) is the ho- 
mogeneous solution of the Langevin equation given in an explicit form in |]TU[. The 
limiting formula is valid in the infinite x limit as derived in [|10|. Both at early and 



late times the factor / ^ 1, deviating from unity only for the time range of order 
1/ujb- 

For discussion of a more general case of finite ^(oo) < we use the local, Ohmic 
approximation, which becomes excellent at late times. A potential that decreases 
fast at infinity as x — ^> oo is assumed; ^ . The acceleration term can then be 
neg lectedif if^l <^rf. This is a slow rolling approximation, and it always holds for 

^cl 



X large enough. The classical equation is then solved as rj /^^ dz{-^)^ = —t , which 
gives the factor / ^ — / {r]p{oo)) . Thus, the tunneling probability decreases 
with time along with the decreasing slope of the potential. Since the tunneling may 
not be terminated within a limited finite time, this result poses a curious question; 
there is a situation in the early universe in which the tunneling is never ended if the 
friction is strong. 

We next turn to the case of V{oo) > Vm, and consider the thermal initial state. 
At low temperatures of T <^ Vh the dominant contribution comes from the en- 
ergy range < E < V^- Contribution from the subbarrier region in integration 
is small and may be ignored. We add two contributions from the left well and 
from the overbarrier region. In the left well region the Wigner function is given by 
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fwiO-i jPi) exp[— tanh(/3a;*/2)(a;*g? +pf/a;*)] . It is more appropriate here to com- 
pute the transition probabihty by ~ / dp fl^^ {x ,p;t) , without the flux p factor. 
The flnal position x is taken around Xq of width Ax. If one makes a rough approxi- 
mation for both the left well and the overbarrier region by harmonic oscillators, the 
relevant classical motion is given by a linear function of initial values, Qi ,Pi, reducing 
the phase space integral to a restricted Gaussian type. 

The averaged tunneling probability over Ax (Gaussian weight assumed) is 
exp (- ^) from the left well region, with A = {^+uj^g^) coth(^) + 2{Ax'^ + 



"-""i" V A J ------------- ...... .^t,—-) V^, I -*y27-— -V 2 

III) , where gi{t)'s are two independent solutions to the linearized equation in the left 
well. The time averaged value over At ^ gives gf ^ and g2 ~ e^^'*/ (2c<j^), 

in the weak coupling. In the late time limit In ~ coth(/3a;*/2)/(2cj*) , and the 
exponent factor becomes of order, exp[— a;4,XQ/(l -|- 2u!^Ax^)] at T ^ a;*, while at 
T > cu, it is exp[-u}lxl/{2T+2u}lAx^)] . 

On the other hand, one may use the WKB form of the wave function in the 
overbarrier region, with the transmission coefficient \T(E)\'^ — Q-'^'^{^h-E)/LOB ^ using 
the inverted harmonic oscillator of the curvature ujb for the subbarrier region. The 
density matrix p^q\q ,q') is then (T » a;* is assumed) 



exp 

Jo 



^{q;E)^*{q'-E), (18) 



where (p{q ; E) cos(/^^ dxp{x ■,E) — j)/ \Jp{q ; E) . The contribution from the over- 
barrier region is then roughly 



/?6,;*(e-2'^'^''/'^«-e-'^^'^)/( {(3-27r/ujB){(^o^ ^ttAx'^ + 7rcoth{(3uJo/2)/uJo)) . In the case 
of overdamping this asymptotic value is approached from below by a slowly varying 
term oc e^^'^o*/'' . Since uj^..Xq ^ Vh and cj* = 0[c<jb], these contributions from the over- 
barrier region are generally larger than those from the left well region. One should 
compare this with the zero temperature result, u^e~'^'"^'^^'^'^ / {uJo^ 2t^Ax^ -|- ti/ujq) 
in order to assess the efi^ect of dissipative medium. 

In summary we gave a real-time formulation of tunneling dynamics in thermal 
medium within the semiclassical framework. 
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Figure caption 



Fig.l 

Two types of tunneling potential. 
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